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Abstract 


Rotating Radio Transients (RRATs) are a relatively new subclass of pulsars that emit detectable radio bursts 
sporadically. We analyzed 10 RRATS observed using the Parkes telescope, with eight of these observed via the 
ultra-wide-bandwidth low-frequency (UWL) receiver. We measured the burst rate and produced integrated profiles 
spanning multiple frequency bands for three RRATs. We also conducted a spectral analysis on both integrated 
pulses and individual pulses of three RRATs. All of their integrated pulses follow a simple power law, consistent 
with the known range of pulsar spectral indices. Their average spectral indices of single pulses are —0.9, — 1.2, and 
— 1.0 respectively, which are within the known range of pulsar spectral indices. Additionally, we find that the 
spreads of single-pulse spectral indices for these RRATs (ranging from —3.5 to +0.5) are narrower compared to 
what has been observed in other RRATs. Notably, the average spectral index and scatter of single pulses are both 
relatively small. For the remaining five RRATs observed at the UWL receiver, we also provide the upper limits on 
fluence and flux density. In addition, we obtain the timing solution of PSR J1709-43. Our analysis shows that PSRs 
J19194-1745, J1709-43, and J1649-4653 are potentially nulling pulsars or weak pulsars with sparse strong pulses. 
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1. Introduction 


Pulsars serve as cosmic lighthouses and can be a powerful 
tool for numerous astronomical studies. They are commonly 
utilized in testing General Relativity (GR; e.g., Kramer et al. 
2006; Antoniadis et al. 2013). Pulsars are frequently timed for 
use in Pulsar Timing Arrays (PTAs) due to their stable spin 
periods (e.g., Hobbs et al. 2010; Manchester et al. 2013) and 
for probing the nature of the interstellar medium (ISM; Han 
et al. 2018). While the majority of known pulsars were 
discovered through periodicity-based searches (Lorimer & 
Kramer 2012), Rotating Radio Transients (RRATs), a subclass 
of canonical pulsars, were detected through single-pulse 
searches instead of the standard Fourier domain searches or 
conventional folding techniques (McLaughlin et al. 2006; 
Keane & McLaughlin 2011). 

RRATS are notable for their sporadic emissions, with only a 
few pulses being detected each hour? (e.g., Keane et al. 2011). 
Another intriguing characteristic of RRATS is their similarity in 
radio emissions to Fast Radio Bursts (FRBs), which are 
extremely bright and manifest as millisecond-duration burst 
events. Single-instance RRAT pulses appear indistinguishable 
from FRB pulses. Dispersion Measure (DM) is used to quantify 
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the time delay caused by the interstellar and intergalactic 
mediums. Generally, FRBs tend to have larger DM values due 
to their extragalactic origin (Bhandari et al. 2018), while 
RRATS tend to exhibit relatively lower DM values (Dong et al. 
2023). Previous studies have attempted to identify potentially 
misclassified FRBs as RRATS, especially when they have only 
been detected through a single-pulse approach (Keane 2016; 
Rane & Loeb 2016). It has been suggested that magnetars, 
another class of neutron stars with high magnetic fields (Kaspi 
& Beloborodov 2017), may be responsible for high-energy 
phenomena such as FRBs (Bochenek et al. 2020; CHIME/FRB 
Collaboration et al. 2020). Consequently, discovering more 
RRATS and magnetars can provide insights into the distribu- 
tion of Galactic pulsars and potentially help characterize the 
populations of RRATs, FRBs and magnetars. 

RRATs are Galactic pulsars characterized by extremely 
variable emission of single pulses. They emit individual pulses 
but then exhibit long periods with no detectable emission. 
However, only several telescopes with high sensitivity are 
thought to be suitable to detect such intermittent emissions and 
obtain sufficient single pulses and periodic emission statistics. 
This inherent challenge in RRAT detection has resulted in 
incomplete knowledge regarding the fraction of RRATS that 
exhibit nulling behavior. Currently, a large number of RRATs 
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lack a precisely measured rotation period, and approximately 
two-thirds have not had their burst rates determined due to a 
lack of follow-up observations. This lack of data significantly 
hinders our ability to properly characterize these sources, as 
highlighted in the study by McKenna et al. (2023). Conse- 
quently, in terms of RRATS, it is necessary to conduct long- 
term monitoring and timing observations. 

The sporadic emission mechanism of RRATSs remains 
unknown, primarily due to the inherent difficulties associated 
with their detection. Despite the discovery of nearly 170 
RRATs as of 2023 (e.g., Tyul'bashev et al. 2018; Good et al. 
2021), the elusive nature of RRATs has posed significant 
challenges in unraveling their emission behavior. Since their 
initial discovery, several models have been put forward to 
explain these sporadic phenomena. Some examples include: a 
normal pulsar exhibiting extreme nulling (e.g., Wang et al. 
2007; Burke-Spolaor & Bailes 2010), fallback of supernova 
material (Li 2006), trapped plasma being released from 
radiation belts (Luo & Melrose 2007) and interference of 
asteroidal or circumpulsar debris (Cordes & Shannon 2008). 
Alternatively, it could be related to mechanisms within the 
pulsar magnetosphere (e.g., Timokhin 2010; Li et al. 2012; 
Melrose & Yuen 2014). Efforts have been devoted to 
understanding this intriguing phenomenon. To establish 
potential connections between RRATs and normal pulsars, 
Zhou et al. (2023) conducted an in-depth analysis of 76 
Galactic RRATSs and their enigmatic emission patterns. The 
study concluded that RRATs predominantly are weaker pulsars 
with a few strong pulses or extreme nulling pulsars. 

In this paper, we present observations of 10 RRATs 
conducted at the Parkes 64 m radio telescope, utilizing the 
ultra-wide-bandwidth low-frequency (UWL) receiver, which 
covers the continuous frequency range of 704—4032 MHz 
(Hobbs et al. 2020). This receiver system provides unprece- 
dented broadband information, including valuable polarization 
data, that plays a pivotal role in advancing our understanding of 
radiation and various phenomena within the intervening ISM, 
such as scintillation and DM variability (e.g., Karako-Argaman 
et al. 2015; Heald et al. 2020). Apart from its large fractional 
bandwidth, it also has a low system temperature, allowing it to 
operate effectively even in the presence of strong mobile phone 
transmissions (Hobbs et al. 2020). In recent years, the UWL 
receiver has played a crucial role in numerous scientific 
projects related to high-precision pulsar timing, the broadband 
analysis of pulsar profiles, and the detection of new pulsars and 
transient sources. Additionally, we analyze the radiation 
characteristics of two RRATs observed using other receivers 
at the Parkes telescope. 


2. Observations 


A majority of pulsar observations from the Parkes Radio 
Telescope archive have been available at the CSIRO pulsar 
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data archive since 1991. Data sets for a large number of 
observations are made publicly available after an 18 month 
embargo’ (Hobbs et al. 2011). The observations from the 
Parkes telescope have been obtained under project codes, 
mainly P1016 - Instant GRRATificiation. The RRATs in our 
sample were observed using the Multibeam, H-OH, 10/50 and 
UWL receivers (e.g., Manchester et al. 2013; Hobbs et al. 
2020). Most of these observations were conducted using the 
UWL receiver system, with 3328 MHz of bandwidth centered 
on 2368 MHz and folded in real time into 1024 phase bins. The 
receiver system provides continuous frequency coverage from 
704 to 4032 MHz and has a low temperature (Hobbs et al. 
2020). PSRs J1649-4653 and J1709-43 were observed by the 
Multibeam Receiver with 256 MHz of bandwidth divided into 
1024 frequency channels and folded in real time into 1024 
phase bins. Additionally, the observed bandwidth is normally 
512 MHz for the H-OH receiver and 1024 MHz for the 10/50 
receiver. All data were recorded in the PSRFITS data format 
(Hotan et al. 2004). In the course of observing, we used five 
backend systems, which included MEDUSA and the Parkes 
Digital FilterBank systems (PDFB1, PDFB2, PDFB3 and 
PDFB4; Hobbs et al. 2006). 

We investigated 10 RRATs in the CSIRO pulsar data archive 
and all observation data have been publicly available as of 
2021. The observing data for the UWL receiver were acquired 
in search mode, while the data for the H-OH, Multibeam and 
10/50 receivers were acquired in fold mode. A summary of the 
observations is given in Table 1, which lists the fundamental 
information of the telescope and the project. The search-mode 
data sets were subsequently folded by DSPSR (van Straten & 
Bailes 2011), at the rotation periods, and the DMs of the 
RRATS were acquired using PSRCAT (Manchester et al. 2005). 
The PAZ tool in PSRCHIVE (Hotan et al. 2004) was used to 
eliminate narrowband and impulsive radio frequency inter- 
ference (RFT). To further mitigate potential RFI contamination, 
we utilized the PFITS_ZAPPROFILE tool in PFITS to process all 
the data. PFITS is a software package designed for reading, 
manipulating and processing PSRFITS format pulsar astron- 
omy data from both search mode and fold mode." Flux density 
calibration was made through the observation of the radio 
galaxy Hydra A (3C 218), where the calibration method can be 
found in Liu & Yu (2020). The calibration data are then applied 
to the RRAT observations using the PAC tool to flatten the 
bandpass and transform the polarization products to Stokes 
parameters. To transform the measured intensities to absolute 
flux densities, we used PAAS to form noise-free standard 
templates from observations and then relied on PSRFLUX to 
obtain the flux density. In order to obtain the spectral 
information of RRATs, we divided the frequency bands into 
several parts and utilized the PULSAR, SPECTRA software to 
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Table 1 
Observing Parameters for 10 RRATs 
PSR Duration N, Np Burst Rate Backend Frequency Bandwidth Receiver Mode 
(hr) (pulses) (pulses) (h^) (MHz) (MHz) 
J0627+16 0.7 1140 ve nee MEDUSA 2368 3328 UWL Search 
J0628--0909 1.0 3001 61 60 MEDUSA 2368 3328 UWL Search 
J1850+15 0.7 1867 e E MEDUSA 2368 3328 UWL Search 
J1909+0641 2.0 9730 116 58 MEDUSA 2368 3328 UWL Search 
J1913+1330 0.9 3459 a ne MEDUSA 2368 3328 UWL Search 
J1919+1745 1.0 1737 835 847 MEDUSA 2368 3328 UWL Search 
J1928+15 0.3 3012 Us Zr MEDUSA 2368 3328 UWL Search 
J1946+24 0.3 257 MEDUSA 2368 3328 UWL Search 
PSR Time Span Sub-int Length Backend Frequency Bandwidth Receiver Mode 
(Years) (s) (MHz) (MHz) 
J1649-4653 10.8 10/20/30/60 PDFB1/2/3/4 732/1369/3094 256/1024 H-OH/MULTI/10/50CM Fold 
J1709-43 7.2 10/20/30/50 PDFB3 /4 845 /1390/3380 64/256/512 H-OH/MULTI/10/50CM Fold 


Note. We have listed the observing parameters for 10 RRATs, but only PSRs J06284-0909, J19094-0641 and J19194-1745 were measured in terms of burst rate. N, 
represents the number of total pulses observed, while N; is the number of burst pulses. The sub-int length represents the sub-integration length. 


find the best-fitting model and produce publication-quality 
plots (Swainston et al. 2022). 


3. Analysis 
3.1. Identification of a Single Pulse 


The first step in our analysis is to discern the pulses 
originating from the RRATs. The intermittent nature of RRAT 
emissions renders Fourier-based techniques or conventional 
folding search algorithms impractical. Following the method of 
Xie et al. (2022), we employ the single-pulse search method to 
specifically target individual pulses surpassing a predefined 
signal-to-noise ratio (S/N) threshold. Generally, a single-pulse 
phase can be segmented into two components: the on-pulse 
region and the off-pulse region. The off-pulse energy is 
commonly characterized by a Gaussian distribution. We set a 
minimum detection threshold of 5o for each pulse, where c 
represents the standard deviation of the off-pulse region. To 
filter out RFI, we exclusively consider pulses whose times of 
arrival (ToAs) derived from the brightest pulse observed during 
an observation fall within a 5% deviation from the expected 
phase. For the purpose of minimizing potential contamination 
from RFI, each pulse was visually inspected, including an 
examination of its frequency-phase plots. The total number of 
detected pulses, the count of filtered pulses and the burst rates 
for all our target RRATs are listed in Table 1. We only detected 
single-pulses in three RRATs (PSRs J19194-1745, J1909 
+0641 and J0628--0909). 

The integrated profiles of all observations for these three 
RRATS are shown in the upper panels of Figure 1, in which the 
on-pulse and off-pulse regions are shown in filled blue and gray 
areas, respectively. In the lower panels of Figure 1, we display 
single-pulse stacks of 15 single pulses for each RRAT, in 


which the burst pulses are signified by red solid lines. It is clear 
to see that all of these three RRATS exhibit sporadic emissions. 
The burst rate of PSR J1919--1745 is 847 h^, which is much 
larger than that of both PSR J0628--0909 with a burst rate of 
60 h ! and PSR J1909--0641 with a burst rate of 58 A. 


3.2. Profile 


A flux calibration procedure was applied to all of the 
observational data. For the search-mode data, we conducted 
observations on each RRAT either once or twice, with each 
observation lasting approximately one hour. To manage 
individual file sizes, the data were automatically divided into 
multiple files. We utilized all observation files to generate a 
cumulative pulse profile, aiming to enhance the S/N. For the 
fold-mode data, our analysis focused only on the integrated 
profile of PSR J1709-43. The sub-integration lengths across 
103 observations varied between 10, 20, 30 or 50. To analyze 
the frequency-dependent characteristics of the integrated pulse 
profile, we utilized single observational data acquired from 
three frequencies: 732, 1369 and 3094 MHz. Each observation 
file contained 30 sub-integrations, and pulse profiles were 
observed across all the sub-integrations. Therefore, the sub- 
integration length across observations did not exert a significant 
influence on the integrated profile. 


3.3. Flux Density Measurements and Spectral Properties 


We utilized the PULSAR SPECTRA software package devel- 
oped by Swainston et al. (2022), which incorporates robust 
statistical techniques to decide on the optimal fiting model and 
calculate the corresponding spectral parameters. The publicly 
available software package presents a comprehensive array of 
tools dedicated to the systematic cataloging of pulsar flux density 
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Figure 1. Integrated profiles and partial pulse sequences of three RRATs. 


and automated spectral fitting. The use of this software facilitates 
the identification of the most appropriate spectral model. 
Notably, the PYTHON-based software implements the five 
spectral models, including the Simple power law, Broken power 
law, Log-parabolic spectrum, Power law with high-frequency 
cut-off and Power law with low-frequency turnover. They are 
adequate for describing the spectra of the vast majority of 
pulsars, with the simple and broken power-law models being the 
most widely employed. Although these models are morpholo- 
gical, the spectral index may potentially be associated with other 
parameters of pulsars. The spectral fitting routine implements the 
method illustrated in Jankowski et al. (2018). To compare the 
five models, Akaike information criterion (AIC) was applied as a 
comparison metric, evaluating how much information the model 
maintains about the data without overfitting. The model leading 
to the lowest AIC is considered the most accurate in terms of 
describing the pulsar's spectra. A more comprehensive descrip- 
tion of how to use the PULSAR. SPECTRA software package can 
be found in Swainston et al. (2022). We present the results of 
spectral fitting from integrated pulses of several RRATs. 
Furthermore, we apply a similar spectral analysis procedure to 
individual pulses with high S/N. 


4. Results 
4.1. Search-mode Data 


The radio emission from pulsars is affected by propagation 
through the ISM in both frequency and time domains, spanning 


durations from seconds to several hours. This phenomenon is 
recognized as diffractive scintillation. We estimate the scatter- 
ing time, 7,, at a reference frequency of 1400 MHz using the 
empirical relationship to DM acquired by Kumamoto et al. 
(2021) 


7, = 1.2 x 10-5DM?2(1.0 + 0.00194DM2). (1) 


For the wide frequency band of the UWL receiver, we 
applied the typical relation 7, ocv ^ (Komesaroff et al. 1972) to 
estimate the scattering time at 950 and 3500 MHz. 

The value of scintillation bandwidth in MHz, Ava, is given 
by Cordes & Rickett (1998) as 


Ay uec (2) 


We estimated the scintillation bandwidth for three RRATs at 
frequencies of 950, 1400 and 3500 MHz, as shown in Table 2. 
The entire observation bandwidth was divided into several sub- 
bands. For PSRs J1919+1745 and J06284-0909, the minimum 
bandwidth of each sub-band is 100 MHz, while for PSR J1909 
+0641, it is 200 MHz. The estimated scintillation bandwidths 
for PSRs J1919+1745, J1909--0641 and J0628+0909 are all 
significantly smaller than the minimum bandwidth of each 
corresponding sub-band. Therefore, the impact of scintillation 
effects on flux density is deemed negligible in this study. 
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Table 2 
Scintillation Bandwidth and Minimum Sub-band Bandwidth of Three RRATs 
PSR DM Amin Avos0 Av 1400 Av3500 
(pe cm?) (MHz) (MHz) (MHz) (MHz) 
J1919+1745 142.3 100 0.001 0.007 0.3 
J1909++-0641 36.7 200 0.3 1.5 58.6 


J0628--0909 88.3 100 0.01 0.05 2.0 


Note. We have listed the observing parameters for PSRs J1919--1745, J1909 
+0641 and J06284-0909. Aui, is the minimum channel bandwidth. Avoso, 
Avj490 and Avasoo are the scintillation bandwidth responses to 950 MHz, 
1400 MHz and 3500 MHz, respectively (Komesaroff et al. 1972). 


4.1.1. PSR J1919+1745 


We have obtained 835 burst pulses, constituting approxi- 
mately 48% of all detected single pulses. In Figure 2, we 
present the integrated profiles of the entire frequency band 
using 62 observations as well as those of nine sub-bands. 

It is evident that all profiles exhibit narrow double-peak 
characteristics. The width of the observed pulses is 0.6 ms, and it 
shows negligible variation with a changing frequency. The peak 
flux densities of the two components of the total integrated profile 
are similar. At lower frequencies (918-1600 MHz), the latter 
component dominates. At middle frequencies (1601—2368 MHZ), 
the flux densities of both peaks are similar. Toward the higher 
frequencies (2369-4032 MHz), the earlier component shows a 
slightly larger peak flux density. Furthermore, the flux densities 
of both peaks exhibit a decreasing trend as the frequency 
increases. Using the method described in Section 3.3, we found 
that the spectra of both the integrated and individual pulses can 
be characterized by a simple power law, which takes the 


following form 
$,- (=) i 8) 
Vo 


where o is the spectral index, vo is the center frequency and c is 
a constant. 

Figure 3(a) shows that the spectral index of the integrated 
pulses is —0.9, which is flatter compared to the average spectral 
index of most pulsars (—1.60 + 0.03, Jankowski et al. 2018). A 
spectral fitting analysis was conducted on a sample consisting 
of 49 individual pulses with high S/Ns, and all the results 
followed the simple power law. The average spectral index of 
49 single pulses is —1.0, ranging from —3.5 to 0.01. An 
example of a single-pulse spectrum is displayed in Figure 3(b). 


4.1.2. PSR J1909+0641 


During the observation period of 2 hr, we identified 116 
burst pulses, corresponding to a burst rate of 58 A !. The 
integrated profiles, encompassing the entire frequency band 
and six sub-bands, exhibit a distinct narrow feature, as depicted 
in Figure 4. The pulse width across various frequency bands is 
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Figure 2. Integrated profiles of PSR J1919--1745 across nine frequency sub- 
bands. 


~0.08 ms, and the peak flux density aligns with the general 
trend of decreasing with increasing frequency. 

We performed spectral analysis on the integrated pulse as 
well as on seven distinct single pulses. It is noteworthy that all 
spectra exhibited a simple power-law distribution. As affirmed 
in Figure 5(a), the spectral index of the integrated pulse is 
—1.2, which is relatively flatter compared to the average 
spectral index of pulsars. Furthermore, the average spectral 
index of individual pulses is nearly — 1.0, spanning from — 1.5 
to 0.05, as presented in Figure 5(b). 


4.1.3. PSR J0628--0909 


We detected 61 individual pulses in 1 hr, corresponding to a 
burst rate of 60 h~'. The pulse width is 0.11 ms. The data are 
split into seven sub-bands. We input the flux densities of the 
seven sub-bands into the PULSAR SPECTRA flux density 
catalog. The spectrum of an integrated pulse can be 
characterized as a simple power law with a spectral index of 
—1.3, which is relatively flatter than the average spectral index 
of normal pulsars. Furthermore, we implement spectral fitting 
on 22 high S/N single pulses using the same approach. The 
vast majority of the single pulse spectra also follows a simple 
power law, as presented in Figures 6(a) and (b), with an 
average spectral index of —1.0 and a range of —1.9 to 0.5. 
Notably, there is an intriguing single pulse which can be 
divided into 11 sub-bands, and its spectrum can be modeled as 
a high-frequency cut-off power law, as shown in Figure 6(c), 
which takes the form 


S, = 13 ( = z) V < Ve (4) 
Vo Ve 
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Figure 3. Spectrum of PSR J19194-1745. 
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Figure 4. Integrated profiles of J1909--0641 across six frequency sub-bands. 


where a is the spectral index, ve is the cut-off frequency and c 
is a constant. This model exhibits a spectrum steepening or 
interruption at high frequencies. A possible explanation 
suggests that radiation originates in the inner (polar) gap, 
where electrons are accelerated in an electric field increasing 
from zero at the star's surface. In this process, electron 
acceleration reaches a maximum and decreases to zero as their 
velocity approaches the speed of light. All emitted power is 
concentrated within the radio frequency band (Kontorovich & 
Flanchik 2013). 


4.1.4. The Other RRATs with search Mode 


For the rest of RRATs observed with the UWL receiver, we 
could not create reliable profiles. Therefore, we are unable to 
carry out subsequent processing and analysis, as we have done 
for the previous three RRATs. 


Limits on the flux density of a single pulse can be described 
by the equation (Tang et al. 2021) 


OS / Nmin Tyys 
GJ AVN, tobs f 


where Tys = 22 K is the system temperature, G = 1.8 K Jy is 
the antenna gain for the UWL receiver of the Parkes telescope, 
tops is the observation time, o is the root-mean-square (rms) 
noise level and Av is the full 3300 MHz bandwidth. As for 


periodic signals, Equation (5) should be multiplied by , na 


(5) 


Sim = 


where 6 is the duty cycle. Due to the absence of a measured 
pulse width for PSR J1850--15 in previous investigations, we 
assume the pulse width to be | ms and a flat spectrum, and our 
non-detection of signal with S/N above 7 puts a fluence 
limitation of ~13 mJy ms. Detailed parameters regarding the 
limitation of flux density and fluence for the five RRATs are 
listed in Table 3. 


4.2. Fold-mode Data 
4.2.1. PSR J1709-43 


We collected observational data spanning 7.2 yr, from 
modified Julian date (MJD) 55537 to 58168. A total of 103 
observations were conducted, all in a fold mode. Pulse profiles 
were detected in 41 observations, constituting ~40% of the 
overall data set. The majority exhibits continuous pulse signals 
across all sub-integrations, while only seven observations 
experienced nulling. Notably, the highest detection rate was 
obtained at 1400 MHz, accounting for ~56% of the observa- 
tions. Three total intensity pulse profiles of PSR J1709-43 at 
732, 1369 and 3094 MHz are in good agreement which exhibit 
a narrow single-peaked structure, as presented in Figure 7. 
These profiles follow the trends of decreasing width and peak 
flux density with increasing frequency. In addition to the 
analysis presented in Section 3, we performed a timing analysis 
for PSR J1709-43, which enabled us to determine its rotation 
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Figure 6. Spectrum of PSR J06284-0909. 


period and the first derivative of the rotation period. The full 
timing solutions are listed in Table 4. 


4.2.2. PSR J1649-4655 


Similar to PSR J1709-43, we collected data for PSR J1649- 
4653 covering 11 yr, from MJD 54220 to 58222, with a total 
observation time of 10 hr. All the observations were carried out 
in the fold mode. Out of these, 65 observations (~51%) 
detected the burst pulse profiles. Figure 8 illustrates the 
variation of peak flux density of J1649-4653 with time. The 
flux density of PSR J1649-4653 is found to be significantly 
weak, with a value of 0.14 + 0.02 mJy for the detectable pulse 
profiles. 


5. Discussion 


We calculated the burst rates for three RRATs, PSRs J1919 
+1745, J1909+0641 and J0628--0909. These rates were then 
compared with results from the Arecibo telescope (Deneva 
et al. 2009) and the Five-hundred-meter Aperture Spherical 


radio Telescope (FAST, Hsu et al. 2023), as detailed in Table 5. 
Notably, J1919+1745 exhibited outstanding behavior with an 
observed burst rate of 847 h~', which significantly exceeds the 
previously reported rate of 320 A^! from the Arecibo telescope. 
Furthermore, the burst pulses constituted ~50% of all the 
observed individual pulses. Consequently, we suggest that PSR 
J1919--1745 may be reclassified from an RRAT to a nulling 
pulsar. For PSR J1909--0641, the observed burst rate was 
58h , relatively lower than the rate of 67 h~' observed at the 
Arecibo telescope. For PSR J06284-0909, the burst rate 
observed in this study was 60 A !. In comparison, Arecibo 
reported a rate of 141 h (s /N 55), and FAST reported a rate 
of 270 h^! (S/N>7). We collected data from 103 observations 
for PSR J1709-43 in fold mode. Among these, pulse profiles 
are detected in 41 observations (~40%). Pulse signals can be 
seen across all sub-integrations for most of them, with only 
seven observations experiencing nulling; 128 observations 
were collected for PSR J1649-4653 in fold mode, and pulse 
emissions were detected in approximately 51% (65 observa- 
tions). Our analysis suggests that the variations in burst rates 
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Table 3 
Summary of the Flux Density and Fluence Limits of the Single Pulses and Periodicity Search of Five RRATs with Parkes UWL Receiver 


Ren et al. 


PSR Name Period Toss c Assumed Width Width Flux Density Limit (70) Fluence Limit (70) 

(s) (s) (rms) Single Pulses (ms) Periodic Signal Single Pulse/Periodic Signal Single Pulse/Periodic Signal 
J0627+16 2.2 2442 30.0 0.03 0.3 21.4 mJy/7.5 ply 0.6 mJy ms/2.3 Jy ms 
J1913+1330 0.9 3168 34.7 0.2 2.0 38.0 mJy/30.2 ply 7.6 mJy ms/60.4 Jy ms 
J1928+15 0.4 1210 12.6 0.5 5.0 21.0 mJy/42.9 py 10.5 mJy ms/214.5 Jy ms 
J1946+24 4.7 1166 21.8 0.4 4.0 10.6 mJy/19.6 py 4.2 mJy ms/78.3 Jy ms 
J1850+15 1.4 2554 32.4 0.1 1.0 29.0 mJy/18.1 py 2.9 mJy ms/18.1 Jy ms 


Note. The pulse widths for the periodic signals of PSRs J06274-16, J1913+1330, J1928+15 and J1946+24 are from Deneva et al. (2009) and McLaughlin et al. 
(2009). As for PSR J1850+15, its pulse width was not found in previous works. Therefore, it is assumed to be | ms. Additionally, it is assumed that the pulse widths 


of all individual pulses are one-tenth of the pulse width of periodic pulses. 


PSR J1709-43 


PSR J1649-4653 
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Figure 7. Integrated profiles of PSR J1709-43 across three frequency sub- Figure 8. Flux density as a function of MJD for PSR J1649-4653. 
bands. 
, Table 4 pulses could be detected. The different burst rates observed in 
Ephemeis fof ESR J170- PSR J0628+0909 may be attributed to the sensitivity of 
Parameter Value different telescopes. The FAST facility demonstrates the 
Pulsar Name PSR J1709-43 highest sensitivity, corresponding to the highest burst rate. 
R.A. (R.A., J2000) 17:09:46 Second, variations in the duration of observations can exert an 
(decl., a A js Wo ^a influence. The observation durations of the FAST and Arecibo 
ererence epoc. n . 
Data Span (MID) 55537-58168 telescopes were comparatively shorter than that of the Parkes 
Dispersion measure (DM) 228 pc cm? telescope. Lastly, it is worth noting that the burst rates of 


Rotation frequency (Fo) 
Frequency derivative (F;) 


1.11499999997 Hz 
—3.0038 x 107^ s7? 


RRATS may exhibit temporal evolution. 
We analyzed the integrated profiles spanning multiple 


i ded Œ a Nera m frequency sub-bands for PSRs J1919+1745, J1909+0641 
stare aed a (B,) in 10? à and J1709-43. A notable trend was discerned: the peak flux 
Characteristic age (7) 0.59 Myr densities exhibited a diminishing trend with increasing 


can be attributed to three primary factors. First, variations in 
telescope sensitivity and observation bandwidth play a vital 
role. If telescopes with higher S/Ns are employed for 
observations, it is possible that a greater number of burst 


frequencies, consistent with the previous results. 
Subsequently, we carried out an analysis and fitting of the 
flux density in various frequency sub-bands for both integrated 
and single pulses. This analysis was facilitated using an 
automated spectral fitting software, PULSAR SPECTRA. We 
present the spectral analyses of the integrated pulses emanating 
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Table 5 
Comparison of Burst Rates Detected with Three Telescopes 

PSR Burst Rate Frequency Duration 

Parkes Arecibo FAST Parkes Arecibo FAST Parkes Arecibo FAST 

(h') a’) a’) (MHz) (MHz) (MHz) (hr) (hr) (hr) 
J1919+1745 847 320 732-4032 1390-1490 EL 1.0 0.11 
J1909--0641 58 67 732-4032 1390-1490 ED 2.0 0.15 Ut 
J0628--0909 60 141 270 732-4032 1390-1490 1000-1500 1.0 0.30 0.48 


Note. The observational data in this study were acquired using the Parkes telescope. The Arecibo data were sourced from Deneva et al. (2009), while the FAST data 


were obtained from Hsu et al. (2023). 


from three RRATs, PSRs J1919+1745, J19094-0641 and 
J06284-0909. All of them consistently follow a simple power 
law. The integrated results are consistent with previous work, 
which reported that 7946 of pulsars have spectra with a simple 
power law (Jankowski et al. 2018). 

However, there is a dearth of literature available for 
comparing single-pulse spectral index with our results. It is 
remarkable that only one single-pulse spectral fitting of PSR 
J0628+0909 exhibits a high-frequency cut-off power-law 
spectrum with a cut-off at 1625 MHz. Complex spectral fitting 
typically requires à substantial number of flux densities at 
different frequencies. However, our current limitations lie in 
acquiring a significant volume of high-quality single-pulse 
spectral flux density data. This may lead to a predominant trend 
in our fitting results, favoring a simple power-law model. 
Jankowski et al. (2018) summarized that there are three 
scenarios for deviations from a simple power-law model in 
pulsars that can be considered: an environmental origin of the 
observed spectral features, absorption processes in the 
magnetosphere and non-generic emission characteristics. The 
calculated means of the single-pulse average spectral indices 
for the aforementioned three RRATs are presented in Table 6. 

Nevertheless, there is a dearth of literature available for 
comparing single-pulse spectral indices with our results. 
Previous studies by Shapiro-Albert et al. (2018) reported that 
the mean spectral indices for PSRs J1819-1458, J1317-5759 
and J1913+1330 are —1.1, —0.6 and —1.2 respectively. The 
range of single-pulse spectral indices they obtained is from —7 
to +4. Xie et al. (2022) reported that the mean spectral index 
for single pulses is —3.2 for PSR J0139--3336 and showed an 
extensive range from —11.85 to +3.83. The frequency ranges 
for these two observations were relatively narrow, spanning 
288 and 500 MHz, respectively. Meyers et al. (2019) reported 
that the single-pulse mean spectral index for PSR J2335-0530 
is —2.2, ranging from —2.8 to —1.5. The observational data 
were obtained from two frequency bands: one centered at 
154.24 MHz with a bandwidth of 30.72 MHz, and the other 
centered at 1396 MHz with a bandwidth of 256 MHz. 
Consequently, the observation covered a broad frequency 
range. Our results exhibit relatively lower scatter compared to 


Table 6 
Spectral indices of Three RRATs 
Integrated 
PSR Pulse Single Pulse 
Standard 

Mean Deviation Range 
J1919--1745 —0.9 —1.0 1.0 —3.5 to 0.2 
J1909-+0641 —1.2 —1.0 0.5 —1.5 

to 0.05 

J0628--0909 —1.0 —1.3 0.7 —2.0 to 0.5 


the results from Shapiro-Albert et al. (2018) and Xie et al. 
(2022), which may be attributed to the broader frequency 
range we measured (3300 MHz) and the longer observation 
duration. 


6. Summary 


We performed an analysis of the emission properties for 10 
RRATS, obtaining burst rates for three RRATS and integrated 
profiles across multiple frequency bands for three RRATs. We 
also conducted an in-depth wideband analysis of the integrated 
pulse and single-pulse spectral characteristics using the 
PULSAR SPECTRA software package. Our findings indicate 
that the vast majority of these characteristics align with a 
simple power-law model, featuring relatively flat spectral 
indices and a comparatively small scatter in single-pulse 
spectral indices. For five additional RRATs observed with the 
UWL receiver, we provide their fluence and flux density upper 
limits. Furthermore, we derive the timing solution for PSR 
J1709-43. Based on the classification of 76 recently discovered 
RRATS by the FAST telescope in the previous work of Zhou 
et al. (2023), we suggest that PSRs J19194-1745, J1709-43 and 
J1649-465 may be nulling pulsars or weak pulsars with sparse 
strong pulses. 
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